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Abstract 

We consider the form of the chiral symmetry breaking piece of the effective 
potential in the linear sigma model. Surprisingly, it allows for a second local min- 
imum at both zero and finite temperature. Even though chiral symmetry is not 
exact, and therefore is not restored in a true phase transition at finite temperature, 
this second minimum can nevertheless mimic many of the effects of a first order 
phase transition. We derive a lower limit on the height of the second minimum 
relative to the global minimum based on cosmological considerations; this limit is 
so weak as to be practically nonexistent. In high energy nuclear collisions, it may 
lead to observable effects in Bose-Einstein interferometry due to domain walls and 
to coherent pion emission. 
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1 Introduction 



The possibility of producing quark-gluon plasma in relativistic heavy ion collisions is 
an exciting one, especially from the point of view of observing the chiral / confinement 
phase transition/crossover as the plasma expands and cools. Formation of domains 
in a chiral phase transition where the chiral field may not be oriented along the 
true vacuum has been a subject of many investigations recently. Formation of 
a large domain with a Disoriented Chiral Condensate (DCC) has been proposed 
by Bjorken, Kowalski and Taylor pj in the context of high multiplicity hadronic 
collisions. It was argued in that, as the chiral field relaxes to the true vacuum in 
such a domain, it may lead to coherent emission of pions Q . A motivation for this 
proposal comes from Centauro events in cosmic ray collisions Q. In the context of 
quark-gluon plasma, Rajagopal and Wilczek proposed [Q] that the nonequilibrium 
dynamics during the phase transition may produce DCC domains. They argued 
that long wavelength pion modes may get amplified leading to emission of coherent 
pions. 

One difficulty in these scenarios is that one typically expects domains which 
are not much bigger than the pion size Q. Several studies have focussed on the 
possibility of getting a larger domain. Gavin, Gocksch and Pisarski have argued || 
that large domains of DCC can arise if the effective masses of mesons are small, 
while Gavin and Miiller propose || the annealing of smaller domains to give a large 
region of DCC. 

In this paper we consider the structure of the effective potential in the chiral 
model and note that more general possibilities exist for the symmetry breaking term 
than considered in these previous investigations. This allows for the existence of 
a second local minimum of the potential, in addition to the true global minimum, 
leading to the formation of domain walls which interpolate between the two minima. 
These walls are unstable unless the two minima are exactly degenerate. We consider 
constraints coming from cosmology on the parameters responsible for the existence 
of such walls. We find that such constraints are extremely weak. We then consider 
phenomenological consequences of the richer structure of the effective potential of 
the model, especially from the point of view of the formation of DCC. We show 
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that in these more general models large domain walls naturally form but eventually 
disappear, leading to emission of pions from shell-like structures. Bose-Einstein 
interferometry should be able to reveal any such shells @. We further argue that 
large regions of DCC may arise naturally in these models, and may be able to 
account for phenomena like Centauro events. 

There is some analogy to the physics of spin glasses ||. A spin glass is char- 
acterized by a phase space which has a complicated landscape of valleys. As the 
temperature is reduced, barriers between valleys become significant, and the re- 
laxation times become very long. The system may even get trapped in a single 
metastable state for the duration of the experiment. One may observe hysteresis 
in strong magnetic fields that varies with the experimental conditions. Real spin 
glasses have an anisotropy. In hydrodynamic models of Heisenberg spin glasses one 
is naturally led to consider a coarse-grained free energy which has similarities to the 
effective potential given in eq. (20) below. 

The paper is organized as follows. In sect. 2 we consider the details of the 
effective potential, both at zero temperature and at finite temperature, and discuss 
various possibilities for the symmetry breaking term. In sect. 3 we consider con- 
straints on the parameters of the model arising from the formation of domain walls 
in the early Universe. In sect. 4 we consider potentially observable effects in high 
energy collisions, namely, pion production from domain walls. In sect. 5 we present 
some conclusions. 



2 The Effective Potential 

The sigma model is thought to represent the long wavelength limit of QCD 
[10, Much has been written about the model, both at zero temperature |T2] , [Hf ] 



and at finite temperature [14]. Despite this, the shape of the effective potential at 
both zero and finite temperature may have a nontrivial structure when the pion mass 
is nonzero, which feature seems to have escaped attention; at least it goes against 
the conventional considerations. In this section we explore the effective potential, 
and in the following sections we shall explore its phenomenological consequences. 
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2.1 The vacuum 

We use the sigma model in its linear representation. The Lagrangian is expressed 
in terms of a scalar field a and the pion field 7r. 

C = \ (d,a) 2 + \ (d^f - ± (a 2 + tt 2 - c 2 /A) 2 - V SB . (1) 

The piece of the Lagrangian which explicitly breaks chiral symmetry is Vsb- In the 
absence of this term, the potential has the shape of the bottom of a wine bottle. 
Chiral symmetry is spontaneously broken in the vacuum, the pion is the massless 
Goldstone boson, and the a meson gets a mass on the order of 1-2 GeV. The axial 
vector current, defined by 

A M = (<9 M cr) 7v - (<9 m tt) a , (2) 

is conserved. It is believed that chiral symmetry is restored by a phase transition, 
which is perhaps of second order, at a critical temperature T c ~ 160 MeV; we come 
back to this point later. 

Any quantitative description of Nature at finite momentum and energy must 



include the vector mesons [15]. Since we are proposing here a qualitatively new 
phenomenon we shall neglect them, as well as strangeness. 

The up and down quarks, while very light, are not massless; therefore, neither 
is the pion. Historically, there have been three ways to add a symmetry breaking 
term to the linear sigma model JO]]: 

1 ; Vsb = ~fn m2 ( 7 , which is linear in the a field. The PCAC relation 

= Umln (3) 

is realized as an operator equation. 

2: Vsb = ^m 2 ™ 2 , which is quadratic in the pion field. The PCAC relation is 
realized only in the weak field limit. 

3: Vsb = constant x iV"iV, which is quadratic in the nucleon field. The divergence 
of the axial vector current is proportional to the axial vector nucleon current. 

In principle one could imagine that low energy pion scattering measurements 
could distinguish between (1) and (2). However, by its very nature that occurs in 
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the weak field limit, so the measurements would have to be very precise. Also, low 
energy pion scattering is influenced by the tails of resonances, such as the p meson 



[17 1, which would tend to mask the effects of the nonlinear symmetry breaking 
terms. We are not aware of any experimental evidence which prefers one symmetry 
breaking term over the other. In this paper we shall not consider the possibility 3. 

If one insists on an effective Lagrangian which is rotationally invariant and 
renormalizable, then the most general symmetry breaking potential can be written 
as 

4 

Other symmetry breaking terms one might think of adding simply amount to a 
redefinition of the 8 parameters A,c, e n ,5 n . Relaxation of the renormalizability 
condition would allow further terms but, as we shall see, there is already sufficient 
freedom to generate interesting physics. In what follows we will set 5\ = 82 = 0, 
mainly for simplicity of presentation. 

With the symmetry breaking potential as given above, the divergence of the 
axial vector current is 

= -V^ B (a)7r, (5) 

where the prime denotes differentiation with respect to a. It is clear that PCAC 
is an operator identity only if Vsb is linear in a. The ground state of this theory 
ought to occur at a(x) = a gs > and 7r(x) = 0. We must immediately impose 
three conditions; that the minimum of the full potential occur at <7 gs , that PCAC 
hold for small fluctuations about <7 gs , and that the pion has its physical mass there. 
Therefore, 

Vsb(<t&) = (6) 

a gs (A^ s - c 2 ) + ^ B (a gs ) = °> ( 7 ) 
Aa g 2 s -c 2 = ml, (8) 

and so <7 gs = /„-. 

One must still ensure that the global minimum is really obtained when tt(x) = 
0. To investigate this problem, let us expand the fields about an arbitrary point as 
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follows. 



a(x) = v cos d + a' (x) , (9) 
tt(x) = v sin6> + 7r'(x). (10) 

The primes denote fluctuations about the given point. The full potential is 

V(v,6) = j (v 2 -c 2 /A) 2 + V SB {v cos9). (11) 

We now consider several limits for the symmetry breaking piece of the potential. 

Suppose that ei = f^m 2 , and all other e's are zero. Minimizing the potential 
with respect to v at fixed 6 gives 

v{0) ~ /tt(1 + A) , (12) 

where A = -2(m n /m a ) 2 sin 2 (#/2), m 2 = 2Xf 2 + m 2 , and 

V min (0) ~ Knin(0) ~ 2/ 2 m 2 sin 2 (9/2) . (13) 

Corrections to the approximate equalities are of relative order A << 1. This is a 
tilted wine bottle bottom with only one minimum. 

Now suppose that £2 = m 2 and all other e's are zero. Then 

v 2 (9) = f 2 [l -2ml sm 2 d /m 2 ] , (14) 

where m 2 = 2Xf 2 , and 

V min (9) - V min (0) « ^/ 2 m 2 sin 2 9 . (15) 

We neglect a correction of relative order (m^lm a ) 2 . This potential has two degen- 
erate minima located at 9 = and ir\ 

In general, one may expect that Vsb allows for two minima, one at 9 = 
and one at 9 = it. If they are not degenerate then by a simple redefinition of the 
fields we may, by convention, choose 9 = to be the global minimum. We shall 
investigate what limits cosmology may place on the existence and depth of the 
second minimum in the next section; we shall find that the constraint is extremely 
weak. It is quite surprising to us that neither terrestrial experiments nor pure 
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theoretical computations in QCD so far tell us anything about a possible second 
minimum. 

Our statement certainly goes against the conventional point of view as ex- 
pounded in ref. [18], for example, which says that the symmetry breaking poten- 
tial should be linear in the fields, and for three quark flavors should follow the 



(3,3*) + (3*, 3) symmetry breaking scheme of Gell-Mann, Oakes and Renner [19]. 
The main argument seems to be that this is the simplest description which gives 
reasonable low energy phenomenology. There was some interest in this issue in the 
1970's in regard to three flavor physics. It was found that the addition of symmetry 
breaking terms bilinear in the scalar fields resulted in low energy phenomenology 
as good as, or better than, linear terms alone [2C[. These bilinear terms could have 
the structure (3,3*) + (3*, 3) or they could have components of some other group 



structure. A more recent study [21| has found a nonzero coefficient of a bilinear 
term in the three flavor nonlinear sigma model. We will not pursue the three flavor 
world in this paper, but it is certainly worth doing. 

One might at first think that e n oc mj (where m q is the up or down quark mass) 
so that ei is much greater than €2 and so on. We think it is quite possible that all 
e n oc m q and therefore of comparable magnitude (when scaled appropriately with 
f n ). The argument is that the sigma model is only a low energy effective model of 
QCD and all possible terms which are allowed should be included. Indeed, even if 
one started originally with only a linear symmetry breaking term e'i, loop corrections 
would generate nonlinear terms. These nonlinear terms would have coefficients equal 
to e'i times some function of A and c. Since A and c are big we expect that all e n /" 
would turn out to be comparable in magnitude. 

For clarity of exposition we shall hereafter restrict our attention to the possi- 
bility that only ei and 62 are nonzero. This is sufficient to parametrize the effective 
potential with the freedom to adjust the tilt of the bottom of the wine bottle as well 
as the depth of the second minimum. We have then at our disposal four parameters 
in the effective Lagrangian: A, c, e\ and e%. These parameters must be restricted so 
as to give the proper pion mass, pion decay constant, a reasonable value for the a 
mass, PCAC in the weak field limit, and the condition that the ground state of the 
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theory occur at a = f n and 7r = 0. We obtain 



ml = A/ 2 -c 2 , (16) 
ml = 2A/ 2 + m 2 -e 2 , (17) 
f n ml = e 1 + e 2 U. (18) 

The numerical values chosen in this paper are m n = 140 MeV and f n = 94.5 MeV. 
Lin and Serot [^2| have emphasized that the a meson in this model is not to be 
identified with the exchange of two correlated pions in the isoscalar - scalar channel 
in the nucleon - nucleon interaction. That exchange is rather broadly distributed in 
mass with a peak around 600 MeV. Good phenomenology for low energy pion and 
nucleon dynamics is obtained if the a meson has a mass greater than about 1 GeV. 
For definiteness, we choose m a = 1 GeV. We shall vary e± between and f^m 2 . 
There is no further freedom given the above constraints. 

2.2 Finite temperature 

To estimate what may happen at finite temperature we will calculate thermal fluc- 
tuations to one loop order and furthermore take the high temperature limit. Quan- 
titatively this cannot be very accurate. The relevant coupling constant is large: 
A (m cr // 7r ) 2 /2 50. Qualitatively the result should be alright; however, see 



To proceed, we expand the fields about an arbitrary point, as in eqs. (9-10), 
and determine the masses of the fluctuations. If v points in the third direction in 
isospin space then the eigenvalues of the mass-squared matrix are 

2 \ 2 2 

mi = Xv — c , 

2 \ 2 2 

m 2 = Xv — c , 



r 2 



mi = 2Xv 2 - c 2 - y - W AM - e 2 Xv 2 cos(2#) + ^ , 



ml = 2At» 2 -c 2 -y + ^X 2 v i -e 2 Xv 2 cos(29) + < ^. (19) 

In the high temperature limit of the one loop approximation one keeps only the 
terms of order T and m 2 T 2 . Ignoring terms which are independent of v and 9 we 



get the simple expression 

V(v,0;T) = c 2 + e 2 cos 2 6>--^U 2 -eit;cosfl. (20) 

In the chiral limit one finds, as is well-known, a second order phase transition 
at the critical temperature T c = ^/2c 2 /A = y/2f w . An analysis by Karsch [24] of 
all available lattice simulations of two-flavor QCD extrapolated to zero quark mass 
is consistent with a second order transition with critical indices the same as the 
0(4) model. So at least qualitatively the model and approximations made here 
make sense. However, as emphasized by Shuryak pq| , the sigma model is supposed 
to represent only the long wavelength modes, and certainly does not include the 
contribution from short wavelength modes. For example, as one approaches T c 
from below, the model does not include the 77, p, u, and the whole tower of mesons 
above them. As one approaches T c from above, the model does not include all the 
degrees of freedom represented by quarks and gluons. The energy density of the 
long wavelength modes represented by the pion and sigma degrees of freedom should 
be thought of as sitting on top of a much larger energy density represented by all 
these other degrees of freedom. 

When the up and down quark masses are nonzero chiral symmetry is not exact. 
It cannot be restored at high temperature. If Vsb is an even function of a the 
Lagrangian still possesses a discrete symmetry which is restored at some critical 
temperature. For example, when only e 2 7^ 0, this temperature is y/2(c 2 + m 2 ) /A. 

It is straightforward to show that the zero temperature effective potential has 
a second, local, minimum at 9 = ir which is separated from the global minimum 
at 9 = by a barrier when the inequality e 2 -\/c 2 /A > e\ > is satisfied. [\/c 2 /A 
is just up to corrections of order m^/m 2 .] As the temperature is increased, 
this minimum develops into a saddle at temperature T\ where the curvature in 
the azimuthal direction changes from positive to negative. The saddle eventually 
disappears at the higher temperature T 2 . Here 
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c 2 fei n 2 

(2 



(21) 
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and 




2/3 



(22) 



In figures 1 through 4 we show the evolution of the effective potential with 
increasing temperature for four sets of e. All of these represent a slice through the 
V — a plane, and are normalized such that V m in = at each temperature. One 
can imagine that a system cooling through T\ could get trapped in the metastable 
minimum. We shall consider such possibilities in the following sections. We call this 
the •proximal chiral phase transition, since it is a consequence of the proximity of 
exact chiral symmetry in parameter (quark mass) space, but it is not a true phase 
transition in the thermodynamic sense. 

3 Cosmological Constraints 

Considerations of phase transitions in the early Universe have been very useful 
in restricting particle theory models. We now ask whether cosmology places any 
constraints on the parameters characterizing the effective potential in eq. (1). As 
mentioned earlier, we will be considering only ei and e<i to be non-zero as this is 
sufficient to capture the qualitative aspects of the effective potential. It is well- 
known that when the effective potential has more than one disconnected minima 
then domain walls are produced in a phase transition. For the potential term in 
eq. (4) this happens if e\ = and £2^0. Two regions of space which correspond 
to the two degenerate vacua a = f n and a = —f n , respectively, (see eq. (14)), will 
be separated by a domain wall where the chiral field smoothly interpolates between 
the two vacua. In the context of the early Universe, stable domain walls are almost 
always disastrous unless the phase transition happens at extremely late times. This 
already suggests that the parameter ei cannot be identically zero. 

By expanding a and 7r, as in eqs. (9-10), we can determine the effective La- 
grangian for 9 from eqs. (1) and (20) to be 



Here v minimizes the effective potential at fixed 9 and T; that is, it traces the 




(23) 
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bottom of the valley of the potential. Since this valley has an almost constant 
radius at fixed T, we can neglect its very weak dependence on 9. Hence v = v(T). 

As we will show in the following, cosmology places a lower limit on ei which 
is very small. Thus, as far as cosmological considerations are concerned, we can 
determine the structure of the domain wall by taking e\ = 0. With this, and 
defining 9' = 29, we get the following equation of motion for 9' , 

a9' + e 2 sm9' = 0. (24) 

This is the familiar sine-Gordan equation which is known to have the domain-wall 



solutions [26] 



9'{z) = 4 tan" 1 [exp (y/e 2 z)] , (25) 

-1/2 

the surface energy density ps of the wall is of the order 

p s ~ e 2 v 2 5 = v 2 ^fe 2 . (26) 

When ei is nonzero then, even though domain walls still form, they are not 
stable any more. [It is simple to check that domain walls always form as long as 
e i < fn^2, which is just the condition that there be a local minimum at 9 = ir.] 
Instability of the wall arises because now one minimum (9 = 0) is energetically 
preferred over the other (9 = ir) so the region corresponding to 9 = ir shrinks while 
the region corresponding to 9 = expands. These unstable domain walls then 
disappear in the course of time as the true minimum spreads throughout space. Of 
course, in the context of the early Universe, there is an upper bound on the life-time 
of such walls if they are not to dominate the energy density of the Universe. As 
the standard theory of nucleosynthesis is in very good agreement with observations, 
one would also like that any such unstable domain walls do not influence it. 

A restriction on the life-time of the unstable walls implies restriction on the 
parameters e\ and e 2 , which we now consider. [Other parameters in eq. (1) are 
not constrained by such considerations from cosmology; the only other topological 
objects in the model are Skyrmions, which are supposed to be nucleons.] Let us 
refer to the difference in the energy densities at 9 = ir and 9 = as Apv- From 
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eq. (23) we have Apy = 2eiv. One expects that the instability of domain walls 
becomes significant when the age of the Universe t* is such that the energy excess 
at the scale t* becomes comparable to the energy of the domain wall on the same 
scale [27|. 



A good way to understand this is to realize that when ei ^ one is actually 
considering a situation similar to a first order phase transition. As domain walls 
keep intersecting each other and pinching off, one can consider closed domain walls 
at any given time. The whole region then looks like bubbles of one phase embedded 
in the other other. The bubbles which have 9 = tt inside will always shrink, while 
the bubbles which have 9 = inside should expand, but only if the size of such 
bubbles is larger than a critical size R c . This critical size can be determined by 
considering the total energy Er of a bubble of radius R, 

Air 

E r = -—R 3 Ap v + AttR 2 P s . (27) 

The critical size R c is determined by requiring that Er be stationary with respect 
to small variations in R. Bubbles larger than R c , which have 9 = inside, will 
expand while those smaller than R c will collapse. We find that 

R c = %*-=v&. (28) 

So far we have neglected the fact that the Universe is expanding. If R c , at 
a given time, comes out larger than the horizon size, at that time, then all the 
bubbles with sizes smaller than the horizon will shrink. Bubbles with sizes equal 
to or greater than the horizon always stretch along with the horizon due to the 
expansion of the Universe, irrespective of which phase is enclosed. Therefore, in 
order that the difference between the global minimum and the metastable minimum 
becomes important, R c must be smaller than the horizon size H (corresponding to 
the age of the Universe at time t*). Using the above equation, this implies that 

H > . (29) 

Using the constraint given by eq. (18) we can write this as 

H> r V f 2 = -*!//*. (30) 

U(m% - e 2 ) ei 
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The horizon size at the time the Universe passed through the chiral/confinement 
phase transition/crossover is about 10 6 cm, which is very large compared to QCD 
scale of 10™ 13 cm. Unless £2 is extremely close to the above inequality will 
be trivially satisfied. [Note that the domain walls do not really have to disappear 
much before the time of nucleosynthesis. Therefore the real constraint is somewhat 
weaker than this.] Since v < f w we find the lower bound to be 

ei > 3 x 1(T 13 MeV 3 . (31) 

As long as e\ is larger than this the domain walls will disappear very quickly and 
will not affect the Universe in any significant way. Because the horizon is very small 
at that time, any density fluctuations generated by collapsing domain walls will also 
get wiped out quickly. 

Clearly, the constraint on ei given by eq. (31) (and the corresponding constraint 
on €2) is extremely weak. Since eq. (1) describes an effective theory anyway, it is 
safe to say that cosmology imposes no practical constraints on the parameters of 
this model. 

4 High Energy Nuclear Collisions 

We now consider chiral symmetry breaking in the context of quark-gluon plasma 
formation in a heavy ion collision and the influence of misaligned chiral condensate 
on pion production. As we mentioned in the introduction, the possibility of coherent 
pion emission from extended domains is very interesting and has been a subject of 
many investigations recently ^, [|, ||] . One of the problems in getting a clean 
signature is that such domains are expected to be very small ||. However, these 
investigations have been restricted to the case when only ei is non-zero. As we 
have discussed earlier, there does not seem any reason to exclude other symmetry 
breaking terms in the potential. The structure of the system drastically changes 
when we consider £2 also non-zero, as exemplified by the presence of domain walls. 
In this section we will consider what sorts of signatures one can expect for these 
more general possibilities for the symmetry breaking terms in the effective potential. 
Let us first consider pion production due to different regions of misaligned 
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condensate. Since previous studies have considered only e± non-zero, we first briefly 
comment on this case. It has a uniformly tilted potential with unique minimum at 
a > 0. It has been suggested earlier j| that, in a rapid phase change, the chiral 
field could roll down to different minima in different regions. In a given region, a 
pion condensate could form if the chiral field points in a direction different from 
the true minimum. As this pion field relaxes to the true minimum it will lead to 
coherent emission of pions. 

One difficulty with getting a clean signal in this scenario is that the typical 
domain one expects is very small, of order 1 fm ||. However, we would like to point 
out that this does not exclude the possibility of the formation of pion condensate in 
a large region. For example, consider two adjacent domains where the chiral field 
points in two different directions, say 0\ and 02- As the chiral field relaxes, one may 
expect that both 9\ and O2 will approach zero. However, this really depends on the 
values of these angles. For example, assume that Q\ = tt + a\ and O2 = tt — 02, 
where both a\ and 0:2 are small. Then in the region where the two domains are in 
contact, the chiral field will have to smoothly interpolate between the two angles, 
and hence somewhere in that region it will point in the direction tt. As the outer 
regions of the domains relax to 6 = 0, the chiral field may start to cover a larger 
region of the order parameter space M. . [We use M. to denote the manifold defined 
by the bottom of the valley traced by the minimum of the effective potential.] This 
means that the evolution of the chiral fields in a collection of domains are not 
totally uncorrelated. Essentially, the chiral field defines a smooth map from the 
region covered by the domains into the manifold M. The image of this map is 
actually a smooth and connected patch in the manifold M. Smoothness of this 
patch is enforced by the condition that the chiral field must smoothly interpolate 
in-between any two adjacent domains. When the chiral field relaxes, this patch 
deforms as its portions slide down to the global minimum. It is quite conceivable 
that this type of evolution of the chiral field leads to a pion condensate pointing 
in a direction which is in some way the average direction defined by the patch. It 
is thus not clear that small individual domains imply no pion condensate in large 
regions. 
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We will not, however, pursue this line of argument since €2 = is unnecessarily 
restrictive. Rather, we will show that a non-zero value of £2 leads to the existence 
of large domains in a very natural manner. 

As discussed in sect. 2, non-zero values of £2 lead to the existence of domain 
walls if the inequality e\ < f n €2 is satisfied. [If not then there is no qualitative 
difference between this case and the case with £2 = 0.] These walls are unstable 
for all non-zero values of e±. Consider the chiral phase change in a region of the 
quark-gluon plasma. As now there are two minima, a local minimum at a = v(ir) 
and the global minimum at a = v(0), the chiral field can relax to either of these 
two values. One therefore expects a domain pattern as shown in figure 5a. For 
pion production, the initial size of these domains is not crucial to our model; this 
clearly distinguishes our case from previous considerations, where the size of the 
initial domain was crucial. In figure 5a we have denoted different domains by the 
angle to which the chiral field relaxes in that region. 

Initially, when the temperature is high, the global minimum is only very slightly 
preferred over the local minimum. At those early times, small walls may collapse 
but large ones may simply be stretched by the expansion of the plasma. This is 
very similar to the situation in the early Universe. As the plasma cools, the energy 
density decreases, eventually reaching a point when the walls become unstable in 
the sense that the 6 = minimum becomes favorable over the 6 = it minimum. If 
the expansion is slow then all the walls with sizes smaller than the critical size [as 
given by eq. (28)] will shrink and disappear. In any case, one is led to a hierarchy of 
sizes of collapsing walls (see figure 5b) which will lead to emission of pions. As the 
walls collapse, the average domain size will increase. Once the typical domain size 
becomes larger than about R c , the instability of walls will become significant. After 
this the regions with 9 = will expand and the regions with 8 = ir will contract. 

The simplest prediction is the formation of walls with size of order R c carrying 
excess energy density. It is important to note that these large walls will form 
irrespective of the size of the initial domains. These walls may expand, converting 
the false minimum into the true one, or they may collapse if they enclose the false 
minimum. One would expect generally that a large wall (comparable to the size of 
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the system) may be left enclosing the true minimum as shown in figure 6. Between 
this wall and the outer boundary of the plasma region the chiral field is in the 
metastable phase. As the 9 = phase expands both from the outside and from the 
inside the two walls will meet. This will lead to a shell of the size of the system 
containing high surface energy density. All the energy contained in the two walls 
will be converted to pions. 

The most important feature of these pions is that they ought to be emitted from 
a shell-like region; studies of Bose-Einstein correlations of pions should be able to 
reveal such a shell structure. Investigations in |Q| could be useful in this context. A 
second feature is the possibility that the pions emitted from such a shell, or from 
the collapsing walls, may be coherent. We now address this possibility. 

Consider a closed wall which bounds a region of the false minimum embedded 
in the true minimum. This wall will collapse and all the energy contained in the 
wall will be converted into pions. The initial structure of the wall is determined by 
the details of the variation of the chiral field from the 9 = region to the 8 = ir 
region. We remind ourselves that we are actually dealing with the minimum of the 
effective potential which is topologically a three-sphere S 3 . Different portions of the 
wall correspond to different trajectories which the chiral field traces from the north 
pole of this S 3 to its south pole. Clearly there is no reason to expect that different 
portions of the wall will all correspond to the same curve on S 3 initially. One can 
then think of the initial distribution of the chiral field in the entire wall as a thick 
strip joining the north and the south poles of S 3 . 

As the wall collapses, one will also expect the thickness of this strip to decrease 
because a thick strip costs gradient energy. Since the pion is very light, the shrinking 
of a strip on S 3 may proceed faster than the wall shrinks, especially for the walls 
which do not shrink initially, either due to the expansion of the plasma, or due 
to being least unstable. If this strip on S 3 shrinks significantly before the wall 
completely collapses, then the entire collapsing shell will correspond to a chiral field 
which, though still interpolating between 9 = and 9 = ir, goes through a unique 
plane. For example, this curve may lie entirely in the tt\ — 112 plane. This might lead 
to a pion condensate, culminating in all of the wall energy being emitted in coherent 
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pions. [Of course, the process of shrinking of the strip on S 3 also will produce pions 
but they will not be coherent.] The same consideration can be applied to the type of 
situation in figure 6 where coherent pions may be emitted from a shell-like structure. 

Let us make a rough estimate of the energy contained in the walls. As an 
example, take €2 = 0.8 m^. Then with eqs. (18) and (28) we get the size of the 
critical bubble to be R c ~ 4.5 m" 1 . With the surface energy density of the wall as 
given in eq. (26) the net energy contained in the wall is about 14 GeV. This is a 
very large energy which can lead to a high multiplicity of pions. In this estimate 
we have only considered the energy of the wall, neglecting the difference in energy 
of the two minima. During the collapse (or expansion) of the wall, false minimum 
energy will be converted into the kinetic energy of the wall which should be included 
to get the net energy emitted in pions. 

We briefly discuss the possibility that our model can also account for Centauro- 
like events |H|. A highly energetic cosmic ray collision may produce a tiny bubble 
of false vacuum such that the bubble wall propogates initially outwards due to 
the initial momentum (or may be due to the initial expansion of the partons). 
Eventually this wall will collapse back. As this wall first expands and then collapses, 
there may be enough time to develop a pion condensate on the wall (due to shrinking 
of the strip on 5 3 in the sense described above), and hence lead to the emission of 
coherent pions. This is especially likely as the initial size of the bubble may be very 
small so the pion field configuration in its wall may be pretty much uniform any 
way. Another possibility is that as the bubble of false vacuum expands due to initial 
wall momentum, a true vacuum bubble of critical size nucleates inside it. This then 
may lead to a large spherical shell containing high surface enrgy density (similar to 
that in figure 6). All this energy may then be emitted in coherent pions. 

If we assume that all of the coherent pions come from a 9 = n bubble, which 
eventually collapses, then the energy emitted in pions E n can be related to the 
radius Rf of this bubble using eq. (27) as 



Rf here represents the radius of the bubble at its largest size. The entire event will 
include any hadrons produced when the bubble was nucleated as well as the hadrons 




(32) 



17 



emitted at the end when the bubble completely collapses. Presumably coherent 
pions will be emitted at the end of the event. If we assume a relation between Rp 
and the duration of the event r, say Rf ~ T a , where a is some parameter, then 
our model predicts a very specific dependence of E n on r (given by eq. (32)). If 
the information of E n and r is experimentally available, then this equation can be 
fitted with data to check our model and hopefully get the parameters a, ei and 62- 

5 Conclusion 

Our main idea is that the effective potential, or free energy, resulting from chiral 
Lagrangian models of QCD may have a second metastable minimum at a chiral 
angle of ir. We illustrated this in the linear sigma model with two quark flavors. 
Our numerical examples were restricted to a symmetry breaking potential which 
had terms linear and quadratic in sigma but, in general, there is no reason to 
think that the other possible terms are ignorable. There does not seem to be any 
fundamental reason why one should not take this situation seriously. Indeed, three 
flavor models with linear plus bilinear terms were investigated briefly in the 1970's. 
The motivation then was to investigate the pattern of symmetry breaking, the 
analytic behavior of observables as the symmetry breaking parameters were sent to 
zero, and to obtain improved phenomenology. It should be kept in mind that the 
parameters in a two flavor model can be renormalized by the presence of a heavier, 
third flavor. Generalizations of our study to three flavors should be done. 

Cosmology places a constraint on the height of a possible second minimum 
relative to the true minimum. This constraint arises from the requirement that 
the energy in domain walls not upset standard calculations of nucleosynthesis. The 
constraint is so weak that it has no practical consequences for high energy particle 
or nuclear experiments. 

High energy nuclear collisions seem to present a remarkable opportunity to 
study the topography of the effective chiral potential at finite temperature. Nonzero 
up and down quark masses spoil the ideal chiral symmetry and smear out the 
probable second order phase transition. This may be a cloud with a silver lining if 
a second metastable minimum exists as it could mimic the effects of a first order 
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phase transition. We have argued that formation and evolution of domains, with 
their attendent domain walls, can plausibly lead to observable consequences. These 
include coherent pion emission and Bose-Einstein interferometry of shell structures. 
Detailed predictions with a specific effective potential require numerical simulations 
as well as the inclusion of vector mesons. 
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Figure Captions 



Figure 1: Temperature dependence of the effective potential for the choice of pa- 
rameters €\ = €2 = 0. The curve shown represents a slice through the V — a plane. 
The potential is rotationally symmetric. There is a second order phase transition 
at T c = y/2f n = 133.6 MeV. 

Figure 2: Temperature dependence of the effective potential for the choice of pa- 
rameters ei = f^rn^ and €2 = 0. The curve shown represents a slice through the 
V — a plane. The bottom of the potential is tilted. There is a saddle point at 9 = tt 
for T < 114.9 MeV. 

Figure 3: Temperature dependence of the effective potential for the choice of pa- 
rameters €i = and €2 = m^. The curve shown represents a slice through the V — a 
plane. There is a second order phase transition restoring the discrete symmetry 
a — ► —a at T c = V2f n = 133.6 MeV. 

Figure 4: Temperature dependence of the effective potential for the choice of param- 
eters ei = 0.25 f-n-m^ and €2 = 0.75 m^. The curve shown represents a slice through 
the V — a plane. The direction 9 = it has a local minimum when T < T\ = 123.2 
MeV, a saddle point when Ti < T < T 2 = 127.0 MeV, and no critical point at all 
when T > T 2 . 

Figure 5: (a) Formation of domains after the phase transition. Domains denoted 
by and tt here represent regions where the chiral field has settled to the true and 
the false minimum, respectively. Solid lines separating different domains show the 
initial structure of domain walls. Outermost solid line denotes the boundary of 
the system, (b) As domain walls join and collapse a hierarchy of domain sizes is 
generated. The solid lines again represent domain walls separating different minima. 

Figure 6: Eventually one may be left with a large shell-like domain of metastable 
matter. As this domain shrinks, one will be left with a large shell (of the size of the 
system) containing large surface density. 
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